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EFFECT CF STANDING TRANSVERSE ACOUSTIC OSCILLATIONS ON FUEL- OXIDANT 
MIXING IN CYLINDRICAL COMBUSTION CHAMBERS 


By William R. Mickelsen 


SUMMARY 

The effect of acoustic oscillations on fuel-oxidant mixing is ana- 
lyzed for acoustic conditions commonly found in screeching and screaming 
combustors. Transverse acoustic fields in cylindrical ducts are described 
in terms of pressure fluctuation, particle velocity, and oscillation 
frequency. The effect of such acoustic fields on the fuel-oxidant mixing 
downstream from point sources of injection is treated theoretically to 
obtain expressions relating the fuel-oxidant ratio to the parameters of 
the stream and the acoustic field. Numerical solutions of these relations 
are made to illustrate the effect of sound-pressure level, oscillation 
frequency, stream velocity, and turbulence level. The results of this 
analysis show that, for moderate screech or scream levels, the fuel- 
oxidant mixing wake is given a transverse oscillatory motion of consider- 
able magnitude and that the fuel-oxidant ratio undergoes large cyclic 
fluctuations which are in phase with the pressure fluctuations of the 
acoustic field. A number of possible mechanisms contributing to screech 
and scream in various combustor configurations are proposed, and methods 
for their control are suggested. 


INTRODUCTION 

Tbe effect of acoustic fields on vapor mixing is of interest to the 
combustion problem in screeching or screa mi ng combustors. If the acoustic 
osc illat ions change the fuel-air or fuel-oxidant ratio, the combustion 
process will be directly affected. To illustrate the effect of acoustic 
osci lla tions, an analysis has been made of vapor mixing in cylindrical 
ducts under conditions usually found in screeching ramjets and after- 
burners and in screaming rockets. 

Acoustic oscillations in cylindrical chambers have three general 
modes : longitudinal, radial, and tangential, as illustrated by the 

following pressure fluctuation diagrams showing pressure fluctuation 
contours at an instant in time: 
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Pressure fluctuation maximums alternate in sign with time. The modes 

shown are the fundamentals; higher modes exist for each type. As pointed 

out in references 1 to 4, acoustic oscillations in ramjet combustors, a 

afterburners , and rockets may he any of these types, depending on a 

number of aerodynamic and combustion factors, 

P 

The pressure fluctuations in each of the modes of oscillation give 
rise to acoustic particle velocities which can affect the aerodynamic 
mixing in a manner Bomewhat similar to that of turbulent velocity fluctua- 
tions. This effect has been analyzed in reference 5, where it is shown 
that the action of a transverse acoustic particle velocity produces a 
sinusoidal-shaped deformation of the concentration pattern downstream 
from a source of heat or mass. A typical deformed mixing wake at an 
instant in time is shown in the following sketch: 



ft til 


Transverse velocity imposed by standing acoustic wave 

11 111 


(b) 
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The analysis of reference 5 shows that, if the deformation is severe, 
the concentration of the heat or mass is changed considerably. In addi- 
tion to the transverse motion of the wake, the shape and magnitude of 
the concentration profile at any axial station is shown to change with 
time. The deformation of the wake is shown to increase with particle 
velocity and, hence, with sound-pressure level. The pressure perturba- 
tions reported in references 2 and 4 are sufficiently large to produce 
transverse particle velocities of the same order as the stream velocity. 
Acoustic oscillations of these magnitudes can strongly affect the con- 
centration level. 

To illustrate the effect of transverse acoustic oscillations on 
fuel-air or fuel-oxidant mixing, the analysis in reference 5 is extended 
herein to a range of acoustic conditions which exist in jet-engine and 
rocket combustors. The acoustic field Is described in terms of pressure 
fluctuation and particle velocity patterns and oscillation frequency. 

An idealized mixing model is used which consists of a point source of 
fuel vapor injected at a constant rate into a flawing, turbulent oxidant 
stream. Mixing equations are derived which define the fuel-oxidant ratio 
at any point in the field and at any time during the acoustic oscillation. 
The effect of tangential acoustic oscillations on the vapor mixing is 
shown to be a composite function of acoustic pressure and frequency, 
stream velocity, and turbulence. Each of these four variables is treated 
separately so that trends in fuel-oxidant ratio may be observed 
graphically. 

In addition to the specific conclusions drawn from the results of 
this analysis, a discussion is included which offers mechanisms which 
might contribute to the driving of the acoustic oscillations and suggests 
methods for their control or elimination. 


PROPERTIES OP THE ACOUSTIC FIELD 

Experiment has indicated that screeching ramjet combustors and 
turbojet afterburners and scream i ng rocket combustion chambers have 
acoustic fields approximated by the closed-end cylinder case. With the 
cylinder closed at both ends, the transverse modes of oscillation are 
preserved throughout the length of the combustor. Departure from this 
uniformity of transverse modes is discussed later in this section. As 
a preliminary to the derivation of the fuel-oxidant mixing equations, 
the acoustic field in closed-end cylinders is described in terms of 
pressure fluctuations (sound pressure), particle velocity, and oscillation 
frequency for a range of gas conditions likely to be encountered in 
combustors . 
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Sound Pressure 

The general solution of the wave equation for the pressure fluctua- 
tion is given in reference 6 (p. 398) and includes all possible modes 
(longitudinal, radial, and tangential) and their harmonics. The class 
of oscillations of interest to the present analysis is the tangential 
modes called the pure standing, pure spinning, and combined standing 
and spinning. The pressure patterns of these three tangential modes are 
shown in sketch (c): 



First pure- First pure- Combination 

standing mode spinning mode spinning-standing 


(c) 

The equation for pressure fluctuation used in reference 4 shows that the 
combination spinning-standing mode can be described by two spinning modes 
of unequal amplitudes rotating in opposite directions. Similarly, the 
pure-standing mode can be represented by two spinning modes of equal 
amplitude rotating in opposite directions. Experimental evidence sug- 
gests that the first spinning-standing mode usually occurs in symmetrical 
combustors, while the first standing mode occurs in combustors having a 
sufficient degree of combustion assymetry. For simplicity, only the 
first pure-standing tangential mode will be. considered. 

As shown in appendix B, the pressure fluctuation field in the first 
pure-standing mode may be described in terms of the root-mean-square 
pressure fluctuation: 



( 1 ) 
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where Vp* is the root-mean-square pressure fluctuation, at the cylindri- 
cal coordinates (r ,6), is "tlie root -mean-square sound pressure at 

the pressure antinodes, and J-j_( ) is the Bessel function of the first 

kind and first order. (All symbols are defined in appendix A.) The 
coordinate system is 



The pressure parameter Vp^yfpj| cos 9 given by equation (l) is shown 
in figure 1 as a function of the fraction of duct radius r/R. 


Particle Velocity 

Equations are derived in appendix B which relate the transverse par- 
ticle velocity components v r and Vq to the sound-pressure level of the 

acoustic field. In order to illustrate the relative magnitude of the par- 
ticle velocity across the duct, equations are derived in appendix B for 

the root-mean-square particle velocity components and 


V? - °* 859 S J 0 (i.841 g) - j 2 £.841 


cos 9 


0.9335 - 
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The root -mean-square particle velocity components expressed by equations 
(2) and (3) are shown in figure 1. 

At the duct axis, the resultant transverse particle velocity is 



0.859 ~ 

r 



( 4 ) 


This resultant 
sound-pressure 


particle velocity is shown in figure 2 for a range 

if, 


levels 


Pq and a range of values of a and 


of 

r. 


As can be seen from figure 2, the particle velocities at the duct axis 
are of the same order of magnitude as the stream velocity for the condi- 
tions commonly encountered in screeching or screaming combustors. 


In summary, the magnitude of the root -mean-square particle velocity 
at the duct center can be determined from figure 2 for particular sound 
field and ambient stream conditions. The variation of root -mean-square 
particle velocity in the duct cross section is illustrated in figure 1 
where particle velocity profiles are shown along the duct radius. 


Oscillation Frequency 


The frequency of acoustic oscillations in the first pure-standing 
mode is given by the expression (ref. 6, p. 398) 


/ 


0.295a 

R 


( 5 ) 


This frequency is shown in figure 3 for a range of duct radii and for 
a range of the speed of sound a. 


Non-uniformity of Transverse Mode in Realistic Combustors 

In realistic combustors the acoustic field may depart from the ideal, 
uniform transverse mode described In the preceding sections. For example, 
data reported in reference 2 show a considerable attenuation and phase 
shift -upstream of the flame in a connected-pipe 6-inch screeching com- 
bustor. A theoretical analysis described in reference 7 predicts an 
attenuation upstream of the combustion zone, but no phase shift. The 
attenuation in both of these cases is due to the change In the speed of 
sound between the burned and unburned gases and to the lack of a suffi- 
ciently strong closed-end condition at the combustor inlet. The 


fi/.q? 
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application of the present analysis to realistic combustors must be 
l imi ted to those cases where the acoustic field sufficiently approximates 
the pure -standing transverse mode. 


TURBULENT MIXING IN A TRANSVERSE ACOUSTIC FIELD 

The point-source model is developed in this section to obtain mixing 
equations which define the fuel-oxidant (or fuel-air) ratio at points 
downstream as a function of time, gas stream velocity, and turbulence 
parameters. The transverse displacement of the concentration wake by 
the particle velocity is defined for vapor Injection at any point in 
the Inlet station of the simulated cylinder combustor. Finally, the 
transverse displacement of the concentration wake is included In the 
mixing equation so that the fuel-oxidant ratio is defined as a function 
of the stream and acoustic field parameters. 


Turbulent Mixing From a Point Source 


The concentration C at any point (X,Y,Z) downstream from a point 
source emitting vapor at the rate Wf into a gas stream of homogeneous 
turbulence is (ref. 8) 


r° 


c = 


Wf 




ay 


W 


(X-IJt) 2 +Y 2 +2 2 

at 


( 6 ) 


It is ass ume d that the vapor attains the stream velocity immediately 
after injection. In terms of fuel-oxidant or fuel-air ratio f/A 
equation (6) becomes 




P o( f / A ) 


w* 


(4jt) 3 / 2 


0 


^75 


(x-ut) 2 +y 2 +z 2 

4cd 


at 


( 7 ) 


where U is the stream velocity, t the residence time, and <x> the 
turbulent spreading coefficient. 


As discussed in reference 8, the turbulent spreading coefficient 
is a function of the residence time t, the turbulence Intensity 


and the Lagrangian correlation coefficient 3 $. In the present 

analysis, the Lagrangian correlation coefficient will be assumed to • 
have the form 



& » e-Vt* 
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where t* « and Sf is the Lagrangian scale of turbulence. With 

i given in reference 8, the 


this assumption and with the equations for 
spreading coefficient becomes 


O) 




& 


se 


1 + e 


( 8 ) 


The spreading coefficient given by equation (8) is shown in figure 
4 as a function of the turbulence intensity and scale. For convenience 
in physical interpretation, ai is plotted against distance downstream 
x ** Ut, but must be treated as a function of time in equations (6) and 
(7). The figure shows that near the injector the spreading coefficient 


depends only on the turbulence intensity 




and residence time. 


For the lower intensities, the spreading coefficient is relatively in- 
dependent of turbulence scale even at distances far downstream, so that 
in approximation the turbulence scale need be estimated only roughly. 


Injection on Duct Axis 


The analysis in reference 5 provides equations describing the motion 
of a mixing wake imposed by a plane transverse sound wave. The analysis 

assumes that the root -mean-square particle velocity V? is constant 
in the region encompassing the deformed wake. The particle velocity 
profiles in figure 1 show that, if the point source is located on the 
duct axis and if the wake is not displaced too far from the duct axis, 
the condition is approximately met so that the displacement y of the 
concentration wake is given by 


7 = ymax sln 


M ■*» b - #) 


(9) 


where 


ymax 


y/2 V? 


*/ 


( 10 ) 


The ratio y/ y mnv is shown in figure 5(a) as a function of the dimen- 
sionless axial distance /x/u and the fractional cycle fx. The trans- 
verse motion of the concentration wake is cyclic and produces definite 
nodes at whole number multiples of the dimensionless distance /x/u. 
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The maximum wake displacement y WHY , which occurs at the antinodal sta- 
tions, is plotted in figure 5(b) for a range of particle velocities and 
oscillation frequencies. 


An expression relating y^^ to the duct radius R and the sound 

pressure in the first pure -standing tangential mode may he obtained 

by combining equations (4) and ( 5 ) with equation (lO) : 


Tmax 

R " 



( 11 ) 


Comparison of equation (ll) with figure 1 shows that, for moderately high 
acoustic amplitudes (say, ^Pr/Po * ^.2), the assumption of constant 
particle velocity is fairly good. 


The magnitude of the pressure fluctuations encountered by the mixing 
wake during its transverse motion about the duct axis will be of interest 

in later discussion. The magnitude of the pressure fluctuation V? 
at the extreme wake displacement y^y can be estimated by evaluating 

equation (l) at r » y mH y along the line 0 « 0, and substituting y raflY 
from equation (ll) : 



( 12 ) 


The pressure fluctuation at y maY expressed by equation (l2) is shown 

in figure 6 as a function of the pressure ratio ^pj|/pQ and of y. 

From figure 6, it is seen that, for injection on the duct axis, the 
mixing wake encounters only small pressure fluctuations when the sound 
field has a moderate amplitude. 


From the preceding discussion and from the derivation in reference 
5, a mixing equation may be written for injection at the duct axis. This 
expression is an extension of equation (7) and accounts for the deforma- 
tion of the mixing wake by inclusion of the transverse displacement of 
the wake center y: 
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P 0 ( f / A ) 

w f 


r\m 





(X-Ut) 2 +(Y-y) 2 +Z 2 

4a> 


&t 


(13) 


where CD is given by equation (8) and y by equations (9) and (l0). 


Injection at Points Away from Duct Axis 


If the point source of injection is placed at a point away from 
the duct axis, the " solid -body " motion expressed by equation (9) no 
longer obtains. The pressure gradient encountered by the wake exerts 
an alternately compressive and expansive motion in the transverse direc- 
tion so that the gaussian fuel-oxidant profile expressed by equation (9) 
will be altered. This effect can be described by deriving the exact 
expression for particle motion along the line 6 *» 0. The particle 
velocity along 0 » 0 is given by equation (B19) which may be rewritten: 


dr 





sin (27t/Tr) 


(14) 


By separating the variables and integrating over the distance y 2 - y^_ 
travelled by a particle during the time interval there results: 


rys/ R 


7l/z 


d .C i M 


0.66 





: 0 (l.84lg) - J 2 (l.34l2) r P 0 


sin (2rt/ir) dx 


(15) 


If = 0 and /r 2 « l/ 2, equation (13) represents the maximum 

displacement of any particular particle in the direction towards the 
wall: 




^y 2 /R 

a(y/R) 

1.32 I 

Jol 

yi/ E 

1.841 g) - J 2 (l.841 | 

j ~ r Po 


(16) 


Equation (16) expresses the maximum displacement of a particle starting 
from any point y 1 at time t » 0. The integral In equation (16) was 

evaluated numerically to obtain a curve such as shown In the following 
sketch: 


4379 
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05 

r- 

to 


f- 



(e) 


A fluid particle initially at y]_ moves to y 2 during the same time 
interval in which a particle initially at y 3 moves to y^. From the 
sketch it can he seen that, in a given time interval, a particle initially 
at y-i_ at one side of the wake moves farther than another particle ini- 
tially at y 3 on the other side of the wake. To estimate the degree of 
departure from the solid-body motion, the ratio (y^ - y 3 )/ (y 2 - y^) is 
shown in figure 7 as a function of the initial position of the wake 

center (y-^ + y 3 )/ 2, the sound pressure yp^, and the initial separation 




has little effect 


of the particles y 3 - y-^ The sound pressure 

in the range shown, while the initial separation has a large effect. 

For small initial particle separation, the solid-body motion is closely 
approximated for point-source positions as much as halfway to the wall. 


The compressive effect of pressure gradient on the fuel-oxidant 
profile across the mixing wake was investigated further by means of a 
simple model. It was assumed that the mixing from a point source had 
produced a given concentration profile centered on the injector axis. 

Then the acoustic field was assumed to act on the wake, pushing it alter- 
nately toward the duct axis and the duct wall. At the initial position, 
the profile was assumed to have a gaussian form such as given by the 
approximate solution of equation (7) given in reference (8): 


x 2 +(y-y ± ) 2 

P 0 C f / A ) 1 

= £**5 e (i7) 

where y i is the lateral distance from the duct axis to the point of 
injection, cdjt was evaluated at X = 1.25 feet, and U was assumed to 
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be 500 feet per second. For a given value of the sound-pressure param- 
eter 1.32 ^^/Vpq; 'fete profiles nearest the duct axis and nearest the 

duct wall were constructed with the aid of equation (16) and sketch (e) . 

In using this graphical method, it was noted that, although the density 
of a fluid particle changes as it passes along the pressure gradient, 
its fuel-oxidant ratio remains constant. The results of this procedure 
are shown in figure 8 for two different sets of the injector positions 

y^/R and the sound-pressure parameter 1.32 -y^/fpQ. The solid lines 

in both cases represent the initial gaussian profile directly downstream 
of the point of injection. All the dashed lines represent the distorted 
profiles when the wake is farthest from the injection axis in both direc- 
tions. The distorted profiles have been transposed to the injection axis 
to facilitate direct comparison. For the sound pressure and injector 
position chosen for figure 8(a), the effect of pressure gradient on the 
fuel-oxidant profile was small. When the sound pressure was increased 
and the injector position moved farther from the duct axis, a substantial 
change occurred in the fuel -oxidant profile as shown in figure 8(b). 

In summary, when Injection is at the duct axis and when the sound- 
pressure level Is not too high, the transverse motion of the mixing wake 
in the first pure-standing mode is nearly a solid-body motion as described 
by equations (9) and (10) . When injection Is at points considerably away 
from the duct axis and when the sound-pressure level is high, the solid- 
body motion no longer describes the motion of the mixing wake with accu- 
racy. In this case, an accurate description of the motion of the mixing 
wake can be obtained by graphical means, but the shape of the fuel- 
oxidant profile can only be approximated. When injection Is at points 
considerably away from the duct axis, a rough approximation of wake mo- 
tion can be obtained by using the solid-body motion (eq. (9)) and by 
evaluating the particle velocity for that equation at the Injector axis 
from equation (4) . 

With vapor injection at points away from the duct axis, pressure 
fluctuations at the Injector will cause cyclic changes In the fuel-oxidant 
ratio. The mixing equation may be corrected for this effect by noting 
that the fuel-oxidant ratio of a fluid particle leaving the Injector is 
proportional to Pq/(pq + p^) where Pq is the ambient mean density and 
Is the density fluctuation at the injector. This density ratio may 
be Inserted into the mixing equation as follows: 


P 0 ( f A) 


w f 


(4*) 5 / 2 

/ o rv o 


_ (X-Ut) 2 +(Y-y) 2 +Z 2 
4co 

e dt (18) 
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where is the pressure fluctuation at the injector as given hy equa- 

tion (B5) . The effect of pressure gradient may he included in equation 
(18) hy using equation (15) to determine y. Equation (18) is not actu- 
ally used in calculating fuel-oxidant ratios in this report; it is in- 
cluded only to illustrate the effect of pressure fluctuations at the 
injector . 


MIXING IN FIRST PURE- STANDING TANGENTIAL MODE 

In this section the mixing equations developed in the preceding sec- 
tion for vapor injection at a point source are applied to acoustic fields 
and gas streams with properties in the range encountered in screaming or 
screeching comhustors. The effects of sound-pressure level, oscillation 
frequency, stream velocity, and turbulence on fuel-oxidant patterns are 
treated for injection at the duct axis in the first pure-standing tangen- 
tial mode. Injection away from the duct axis is discussed qualitatively. 
The first pure -standing tangential mode is assumed to exist throughout 
the length of the simulated combustor, and the gas stream properties are 
assumed constant throughout the combustor for any particular set of 
conditions . 


Effect of Amplitude and Frequency of Acoustic Oscillations 

To investigate the effect of sound-pressure level and frequency, 
fuel-oxidant ratio profiles were calculated from equation (13) by numeri- 
cal integration with a high-speed digital computer using Simpson's rule. 
The profiles were calculated fen- a fixed axial station X = 1.25 feet 
and for two instants in time /r = 0.5 and /r = 0.75.. The turbulent 

stream properties were fixed at U = 500 feet per second, Vl/u = 0.05, 

and = 0.1 foot, so that = 1.603*10"3 square foot-. The profiles 

were calculated for a range of particle velocities and oscillation fre- 
quencies. The frequencies were chosen so that the axial station X = 1.25 
feet was an antinodal station in each case. (Antinodal stations are 
illustrated in fig. 5(a) and are at /x/U =0.5, 1.5, . . ..) Antinodal 
stations were chosen for the calculations because the wake attains its 
maximum lateral displacement there. The group of fuel-oxidant ratio pro- 
files calculated for /t = 0.5 are shown in figure 9(a) and have t he 
following geometrical relation to the mixing wake: 
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2200 

As shown "by sketch (f ) , Increases in frequency reduce the lateral dis- 
placement and internodal distances of the wake. 

In figure 9(b), the wake motion has progressed one-quarter cycle to 
S'* — 0.75, and the calculated fuel-oxidant profile lies in a plane normal 
to the wake center: 
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Hie profiles were calculated for planes normal to the wake center to 
facilitate comparison with the profiles in the absence of screech (solid 
lines in fig. 9(b)). For a frequency of 2200 cycles per second, the rise 

in fuel-oxidant ratio to the left of center in the curve for = 600 

feet per second is due to the preceding wake wrinkle. 

As the particle velocity is increased, the peak fuel-oxidant ratio 
decreases in all cases shown in figure 9. In addition to this general 
decrease in fuel-oxidant ratio, a cyclic change in fuel-oxidant ratio is 
also shown. (The profiles in fig. 9(a) will be repeated at /t = 0 and 
those in fig. 9(b) will be repeated at /x = 0.25.) As shown in refer- 
ence 5, the fuel-oxidant ratio will be highest in portions of the wake 
having zero slope and lowest in portions having maximum slope. This in- 
dicates that the fuel-oxidant ratios shown for /t = 0.5 (fig. 9(a)) will 
form the upper limit of the cyclic variation, and the fuel-oxidant ratios 
in figure 9(b) will form the lower limit. By comparison of figures 9(a) 
and (b) it can be seen that the magnitude of the cyclic variation in 
fuel-oxidant ratio is greatest for the lower oscillation frequency and 
increases with particle velocity. 

Consideration of figures 5 and 9 shows that a time and position 
relation exists between the fuel-oxidant and pressure fluctuation. When 
the wake is at its maximum lateral displacement, the fuel-oxidant ratio 
at the antinodal stations is at a maximum, and the pressure is also at 
a maximum. This relation is expressed by equation (B26), derived in 
appendix B, which defines the pressure change experienced by the wake 
center during its transverse motion. The pressure fluctuation calculated 
from equation (B26) is shown in figure 10 In comparison to fuel-oxidant 
fluctuations obtained from figure 9. This comparison is made for two 
axial stations: /x/U = 0.5, 1.5, . . . and /x/U = 0.75, 1.75, .... 

At the antinodal stations /X/U =0.5, 1.5, . . ., the pressure and fuel- 
oxidant fluctuations are exactly in phase as shown in figure 10 (a) . At 
stations away from the antinodal, such as /x/U = 0.75, 1.75, . . . 

(fig. 10(b)), the pressure and fuel-oxidant fluctuations are no longer 
in phase, but the magnitude of the pressure fluctuation is considerably 
reduced. 


Effect of Stream Velocity 

The effect of stream velocity on the mixing pattern is shown in 
figure 11 in which fuel-oxidant profiles are shown for the axial station 
X = 1.25 feet for two Instants In time /t = 0.5 and /t = 0.75. Hie 
profiles were calculated for a fixed oscillation frequency of 1000 cycles 
per second and a spreading coefficient of 1.603*10“ 3 square foot. The 
geometry of the concentration wake and fuel-oxidant profile for figure 
ll(a) Is: 
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(*) 


In figure 11(b) the wake motion has progressed one-quarter cycle, and 
the calculated fuel-oxidant profile lies in a plane normal to the wake 
center : 



Center of 

concentration wake 
Stream velocity, ft/sec 
278 
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Comparison of figures ll(a) and (b) shows that stream velocity has 
an effect on fuel-oxidant ratio similar to that of oscillation frequency 
as discussed in the preceding section. The general decrease in fuel- 
oxidant ratio with stream velocity is more marked than that with, oscilla- 
tion frequency. A cyclic fluctuation in fuel-oxidant ratio is again noted 
and is most pronounced for the lower stream velocity. 


Effect of Turbulence and Residence Time 

The effect of turbulence parameters is illustrated in figure 12 
where fuel-oxidant profiles are shown for two sets of turbulence condi- 
tions. The profiles were calculated for a stream velocity of 500 feet 
per second, an oscillation frequency of 1000 cycles per second, a time of 
/t = 0, and an axial station X = 1.25 feet. The peak fuel-oxidant ratio 
has the same percentage decrease as the particle velocity increases in both 
cases . From this observation, it is concluded that for the range of con- 
ditions investigated the turbulence does not affect the mixing in any way 
other than normal turbulent mixing. 

Since the spreading coefficient is also a function of residence time, 
it also may be concluded that, for the range of conditions investigated, 
the residence time, or distance downstream, does not affect the mixing in 
any but the usual, manner. 


Injection at Points Away from Duct Axis 

As discussed in the earlier section titled Injection at Points Away 
from Duct Axis, injection at points away from the pressure node introduces 
a pressure fluctuation effect as expressed by equation (18) . A qualitative 
picture of this effect is illustrated in figure 13 where the fuel-oxidant 
wake is shown for successive intervals of time along with the pressure 
fluctuation carve. The width of the wake represents the magnitude of the 
peak fuel-oxidant ratio. When the pressure fluctuation is negative, the 
fuel-oxidant ratio at the injector is high, but is low when the pressure 
fluctuation is positive. The figure shows that, when the richest portion 
of the wake has reached its maximum excursion towards the duct wall, the 
pressure fluctuation has reached its maximum positive value. A second, 
lesser maximum of fuel-oxidant ratio reaches its maximum excursion towards 
the duct center when the pressure is at its other maximum positive value. 

In this way, the pressure fluctuations at the injector accentuate the 
fuel-oxidant ratio fluctuations previously discussed in the section 
Effect of Amplitude and Frequency of Acoustic Oscillations. 


POSSIBLE MECHANI24S CONTRIBUTING TO SCREECH AND SCREAM 

The complexity of the aerodynamic and chemical processes which occur 
In a combustor has prevented the formulation of a straightforward, all 
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Inclusive description of the screech and scream driving mechanisms. The 
role of pressure fluctuations in driving screech is discussed in refer- 
ence 2 where it is suggested that a periodic change in pressure is ex- 
pected to cause an even larger periodic change in heat-release rate. 
Consideration of the equation for chemical -reaction rate indicates that 
fuel-oxidant fluctuations would also cause fluctuations in heat-release 
rate. For example, consider the second-order reaction rate: 


Reaction rate "C^(f/A)e 


-e/b c t 


(19) 


where C A is the initial concentration of oxidant, E the activation 
energy, R c the gas constant, and T the flame temperature. Hie flame 
temperature is a strong function of the fuel-oxidant ratio so that the 
heat -release curve has the following form: 


Heat- 

release 

rate 



(j) 

The pressure and fuel-oxidant ratio driving mechanisms are only two of a 
number of possible screech and scream driving mechanisms. Other mecha- 
nisms are described in references 9 and 10. 

The results of the present analysis have shown that the fuel-oxidant 
mixing wake downstream of the injectors is affected by transverse acoustic 
oscillations in two principal ways: 

(1) The fuel-oxidant mixing wake is given an oscillatory transverse 
motion which is quite large for moderate screech or scream intensities. 

(2) During its transverse motion at antinodal stations, the mixing 
wake undergoes periodic changes in fuel-oxidant ratio which are exactly 
in phase with the pressure fluctuations. 

With fuel injection at the duct axis, the fuel-oxidant fluctuations 
caused by the transverse acoustic oscillations may be more important to 
the screech or scream driving mechanism than the fluctuations in pressure. 
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As shown in figure 6, the pressure fluctuations encountered hy the mixing 
wake during its transverse motion are quite small with injection at the 
duct axis, "but as shown in figure 10, the fuel-oxidant fluctuations may 
he quite large for this case. These observations and results suggest 
that the actions of the mixing wake may he important to the driving mech- 
anism of screech and scream. To pursue this notion further, possible 
screech and scream mechanisms are discussed in detail in this section, 
and methods for their control are suggested. 

In the case of ramjet combustors and afterburners where injectors 
are positioned upstream of flameholders, the axial distance between the 
injector and flameholder stations may be of critical importance to the 
screech mechanism. As shown in figure 5 the mixing wake may be periodi- 
cally displaced a considerable lateral distance at antinodal axial dis- 
tances from the point of injection. If the distance from the injector 
to the flame front is at or near a wake antinodal distance, then the fuel- 
air mixture might pass alternately back and forth across the flameholder 
and cause significant periodic changes in heat release In the two opposing 
flame fronts. The prevention or control of this type of screech mecha- 
nism may lie in the distance from the injector to the flame front. By 
changing this distance by a fraction of the internodal distance, the time 
and position phase relation between the fuel-oxidant ratio fluctuations 
and pressure fluctuations might be altered to attain a suppressing effect 
on the screech driving mechanism. 

The possibility of contributions to screech by periodic fuel-air 
ratio fluctuations is discussed in this paragraph. Unpublished MCA data 
on full-scale combustor operation, shown in figure 14, show that screech 
occurred at lean over -all fuel-air ratios . In the absence of screech, 
the peaks in the fuel-air ratio profile at the flame front may be expected 
to be at stoichiometric or lower for such lean over-all fuel-air ratios: 
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In the section MIXING IN FIRST PURE -STANDING TANGENTIAL MODE, it was shown 
that transverse acoustic oscillations cause both a general decrease in 
peak fuel-air ratio at any station downstream from the injection station 
and a periodic fluctuation in peak fuel-air _ ratio at antinodal stations, 
which is in phase with the pressure fluctuations. This situation is 
illustrated in sketch (l): 


U 


Spray bar 



Stoichiometric 


f/A 


0 f/A 


Heat- 

release 

rate 



CD 

The fuel-air ratio profile shown by the solid line represents the upper 
limit of the fuel-air ratio fluctuation, and the dotted line represents 
the lower limit. The fuel-air ratio profile shown is assumed to be at 
the flame front, which is at one of the antinodal distances /x/u « 0.5, 

1.5, . . . downstream from the spray bar. Since the mean fuel-air ratio 
profile is lower than Btoichiometric, a fluctuation towards higher fuel- 
air ratio will result in an increased heat-release rate as shown by the 
reaction rate curve above. This increase in heat release due to a rise 
in fuel-air ratio is directly in phase with the increase in heat release > 

due to pressure rise, so that the screech is driven to even higher 
amplitude. If the mean fuel-air ratio were higher than stoichiometric, 
then a rise in fuel-air ratio would tend to decrease the heat-release 
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rate and thus exert a suppressing effect on the screech driving mechanism. 
To obtain a mean fuel-air ratio higher than stoichiometric, the number of 
injection points could be reduced with a higher fuel flow rate through 
each: 



(m) 


An alternate method could be to move the original spray bar 
downstream: 



(*) 
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Another method of suppressing the effect of fuel-air ratio fluctua- 
tions is control of the fuel-flov rate. The stretching of the mixing 
wake by the acoustic oscillations could be compensated by mechanically 
varying the fuel flow rate in proper time -phase relation with the pres- 
sure fluctuations. Overcompensation of the stretching effect would cause 
a phase inversion of the fuel-air ratio and pressure relation, thus 
tending to suppress further the screech driving mechanism. 

The effect of offcenter injection on fuel-oxidant ratio fluctuations 
has been discussed in the second section titled Injection at Points Away 
from Duct Axis, where it was shown that the pressure fluctuation at the in- 
jector produces fuel-oxidant ratio fluctuations in phase with those caused 
by the deformation of the mixing wake by particle velocity. It should 
be noted here that the transverse particle velocity decreases as the 
injection axis is moved away from the duct axis in the first pure-standing 
tangential mode. Since the particle velocity decreases, the magnitude of 
the fuel-oxidant ratio f luctuations will also decrease as shown previously. 
From these considerations it appears that, as the fuel injector is moved 
away from the pressure node, the changes in magnitude of the pressure 
fluctuation and particle velocity act in a manner to compensate each 
other, thus diminishing the effect of change in injector position on the 
fuel-oxidant ratio fluctuations. If this compensation is exact, then 
the fuel-oxidant ratio fluctuations will be of the same magnitude as 
obtained with injection at the pressure node. However, the mixing wake 
will encounter more severe pressure fluctuations during its transverse 
motion so that the screech driving mechanism might be strengthened by 
positioning the injector away from the duct axis. 


CONCLUSIONS 

In summary, equations have been derived for the calculation of fuel- 
oxidant ratio profiles downstream from injectors simulating point sources 
in combustors with transverse acoustic oscillations. 

For combustors having the first pure-standing tangential mode of 
acoustic oscillation and injection occurring on the combustor axis, the 
following conclusions are made: 

1. Acoustic oscillations of the magnitude present in screeching or 
screaming combustors can cause considerable transverse motion and distor- 
tion of the fuel-oxidant mixing wake. 

2. In general, the peak fuel-oxidant ratio is lowered by the trans- 
verse acoustic oscillations. For example, when the peak-to-peak amplitude 
of the pressure fluctuation is equal to the mean c hamb er pressure, the 
peak fuel-oxidant ratio may be reduced as much as 60 percent in a gas 
stream at 4000° R flowing at 300 feet per second in a 2 -foot -diameter 
duct. 
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3. The decrease in peak fuel-oxidant ratio becomes more marked as 
the amplitude of the acoustic oscillation increases. 

4. The decrease in peak fuel-oxidant ratio becomes more marked as 
the stream velocity decreases. 

5. At any particular combustor station, the peak fuel-oxidant ratio 
has a cyclic variation with time. The phase relation between the peak 
fuel-oxidant ratio and pressure fluctuations is a function of the posi- 
tion of the station, the stream velocity, and the oscillation frequency. 
For the cases of low oscillation frequency and low stream velocity or 
both, the portions of the mixing wake which are displaced farthest from 
the duct axis have a cyclic variation in fuel-oxidant ratio which is 
exactly in phase with the pressure fluctuations (i.e., maximums of fuel- 
oxidant ratio and acoustic pressure occur at the same time). 

For combustors where injection occurs at other than the pressure 
node of the first pure tangential mode of acoustic oscillation, the 
following conclusion is made: 

1. With constant fuel injection rate, the pressure fluctuations at 
the point of injection accentuate the fuel-oxidant fluctuations at 
stations downstream. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, March 18, 1957 
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APPENDIX A 
SYMBOIfi 

constant 

speed of sound, ft/sec 

concentration of diffusing vapor, Ib/cu ft 
initial concentration of oxidant 
activation energy 
Naperian "base 

frequency, or frequency of combined transverse and axial nodes, cps 

frequency of combined transverse modes, cps 

frequency of axial mode, cps 

fuel-air or fuel-oxidant ratio by weight 

Bessel function, first kind, order n 

length of cylindrical duct, ft 

Lagrangian scale of turbulence, ft 

wave number for radial mode 

wave number for tangential mode 

wave number for axial mode 

sound-pressure fluctuation, lb/sq in. 

pressure fluctuation at injector, lb/sq in. 

ambient mean static pressure, lb/sq in. 

root mean square of sound-pressure fluctuation, lb/sq in. 


c f r**. 
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R 

R 

c 

9t 

r 

T 

t 

t* 

XT 

v 

v c 

v r 

^0 


root mean square of sound-pressure fluctuation at pressure anti- 
node, Ib/sq in. 

radius of cylindrical duct, ft 
gas constant 

Lagrangian turbulence correlation coefficient 
radial coordinate in cylindrical coordinate system, ft 
flame temperature 

residence time in fluid stream, sec 

characteristic time of turbulence, sec 

stream velocity, ft/sec 

resultant particle velocity, ft/sec 

resultant particle velocity along duct axis, ft/sec 

radial particle velocity, ft/sec 

tangential particle velocity, ft/sec 



root mean square of particle velocity. 


ft/sec 



X,Y,Z 


x 


Ymax 

c ^n,n 


Intensity of turbulence, ft/sec 

fuel flow rate from point source, Ib/sec 

position of particular point in rectangular coordinates, ft 
axial coordinate in cylindrical coordinate system, ft 
lateral displacement of mixing wake center, ft 
lateral distance from duct axis to point of injection, ft 
maximum lateral displacement of mixing wake center, ft 
roots of the equation [dJ m (rta)/d£t] » 0 
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r 

0 

Pi 

p o 

T 

CD 

CDjj 


ratio of specific heats 

angle of cylindrical coordinates , radians 
density fluctuation at point source , Ih/cu ft 
mean ambient density, lb/cu ft 
time scale of acoustic oscillations, sec 
turbulent spreading coefficient, sq ft 

turbulent spreading coefficient at axial station X, sq ft 


4379 



CS-4 tack 4379 

( I 


NACA TN 3983 


27 


APPENDIX B 


DERIVATION OF EQUATIONS 
Pressure Fluctuation 


Hxe general solution of the wave equation for the pressure fluctua- 
tion p is given in reference 6 (p. 398) : 

P - - PO^n ( — a^ ~ ^ cos ( n0 ) cos (~ ~ a ~) e" 12 ^ (Bl) 

where the minus sign is used for convenience in later analysis and where 
the cylindrical coordinate system (x,r,0) is used. The various frequen- 
cies are 


J „ a 

' m^n 


a 


m,n 

2R 


4 



(B2) 


where the values of ct ta ^ n are given in reference 6 (p. 399). Hie my n, 
and n x are wave numbers for the radial, tangential, and axial modes, 
respectively. 


The expression for the pressure fluctuation in the first pure- 
standing tangential mode is 


0,1 \ -i2rt/ n 

- PqAJ-^ l — vj cos 6 e 


cos 6 e (B3) 

' - \ ” r 

which may he written as 

p => - PqAJ-j^ rj cos 6 [cos (Zk/t) - i sin ] (B4) 


Using the real part 


- ^ 1 ^ r ) 


cos 9 cos (2rt/r) 


(B5) 


at the pressure antinode, where 0 b 0 and r «= R, and noting that 
/ =s 0.293 a/R give 


a AJ-l ( 0.5861 jt) 


V2> P B 

PO 


(B6) 
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so that 


A - 2.43 



Substituting equation (B7) into (B5) gives 

p = -2.43 ^p| ^ ^1.841 cos (2it/r) cos 0 

from which 



(B7) 


(B8) 


(B9) 


Tangential Particle Velocity 

For the magnitudes of pressure f luctuations considered herein, the 
particle velocity is related to the pressure fluctuation by the acoustic 
relation given in reference 6 (p. 295): 

v « ; - g ■ - grad (p) 
i2jt/rp 

The tangential component Vq is obtained from equation (BIO) : 



g 


(BH) 


** i2^/rp o9 

Substituting equation (B7) into equation (B4) results in 

p = -2.43 -^p^ Jp ^1.841 cos 0 [cos (2«/t) - i sin (2jt/t)] (B12) 

By noting that / =* 0.293 a/R and g(ap) = a/(yp 0 ) the following equa- 
tion is obtained: 


Vq w 1. 


.32 - J, (l.841 sin 0 [-i cos (2 */t) - sin (2 «/t)] 

r p o r \ r-o 


(B13) 
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Using the real part results in 


1 a t P R 
6 » -1.52 r Po 


I J x jl.841 sin 0 sin (2 ji/t) (B14) 


In terms of root -mean-square particle velocity, equation (B14) becomes 

sin 0 (B15) 

In the limit at r =» 0 and along 0 =* jc/ 2, equation (B15) becomes 

(B16) 


“ 0.9335 S | Ji(l.841 g) 




0 

Radial Particle Velocity 

From equation (BIO) the radial particle velocity component is 

-S — ^ (B17) 


T r 


±2rtfp 


Noting that /=» 0.293 a/R and g/ap » a/rp 0 ani 116 in S equation (B12) 
result in 



- |j 0 ^..841 - J 2 ^L.841 cos 0 [-i cos (2 */t) - 

(B18) 


sin (2it/ar) ] 

Using the real part gives 

it 


V r « 1 y - 1 -- [j 0 841 |j - J 2 ^1.841 cos 0 sin (Z *A) (B19) 

In terms of the root -mean-square particle velocity the equation is 

' 2 r / \ / \n 

cos 0 (B20) 


4 m 0.859 S j 0 (l.841 g) - Jg (l.841 g) 
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and at r *» 0 along 9 » 0 It "becomes 



(B21) 


Fressure Fluctuations Experienced "by Mixing Wake with 
Injection on Duct Axis 

To augment the discussion of results, it is of interest to describe 
the pressure fluctuations encountered "by the mixing wake during its 
transverse motion. The pressure fluctuation experienced by the wake 
center can be expressed by using the wake displacement y in place of 
the radius r in equation (B8): 


p » -2.43 |l. 841 cos (2 k/t) 


(B22) 


By noting that the particle velocity near the duct axis is approximately 
given by equation (4), equation (9) may be written: 


y , q-,859 >/2a JYgB sj[n 

p 0 


(^) sin (zx/t - ^ (B23) 


Substituting equation (B23) into (B22) and noting that /» 0.293 a/R 
give 


5.43 ^ J ± 


2.43 


sin 


(#) ■*> h - $ 


cos (2 ic/t) 

(B24) 


Expanding the Bessel function in its infinite series and noting t ha t its 
argument 1.841 y/B is small (see eq, (ll) ) so that the higher-order 
terms may be neglected give 


1.215 V3 (tt) =1“ ( 2 *A - 


cos (2rtfx) 


(B25) 


c t 
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Normalizing gives the following equation: 


rv qP 

PR 


-2.95 sin 



(B26) 


At antinodal stations, /x/u * 0.5, 1.5, . . . and equation (B26) 'becomes 


TPqP 

Pr 


2.95 cos 2 (2 ji/t) 


(B27) 
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Figure 2. - Relation between particle velocity at duct axis and sound pressure 
in first pure-standing tangential mode. 
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Figure 4. - Turbulent spreading coefficient for range of 
turbulence intensities and scales. 
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Dimensionless transverse position, y/'y. 
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Dimensionless axial position, /X/U 

(a) Position of concentration wake center at various times through one cycle. 

Figure 5. - Motion of concentration wake with point source on duct axis in first pure- 

standing tangential mode. -J 
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Maximum displacement of wake center, y mHV , 
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Frequency, /, cpa 

(b) Maximum mixing wake displacement. 

Figure 5. - Concluded. Motion of concentration wake with point source on duct axis in first 
pure-standing tangential mode. 
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(a) Fraction of cycle, 0.5. 


Figure 9. - Effect of oscillation frequency on mixing in first pure-standing 
tangential mode. Stream velocity, 500 feet per second; turbulence intensity, 
5 percent; turbulence scale, 0.1 foot; axial station, 1.25 feet. 
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(b) Fraction of cycle, 0.75. 


Figure 9. - Concluded. Effect of oscillation frequency on 
miming In first pure-standing tangential mode. Stream Teloc- 
ity, 500 feet per second; turbulence intensity, 5 percent; 
turbulence scale, 0.1 foot; amial station, 1.25 feet. 
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Fraction of cycle, /c 

(b) At stations fXjXS «= 0.75, 1.75, .... 

Figure 10. - Time-phase relation between pressure and fuel- 
oxidant fluctuations on the wake center during its trans- 
verse motion. 
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Lateral distance from duct centerline, Y, ft 
(a) Fraction of cycle, 0.5. 

Figure 11. - Effect of stream velocity on mixing in first pure-standing 
tangential mode. Oscillation frequency, 1000 cycles per second; 
turbulence intensity, 5 percent; turbulence scale, 0.1 foot; axial 
station, 1.25 feet. 
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(a) ^v^/U = 0.05; 0.10 foot; ax^ = 1.603*10 -3 square foot. 

*° 6 1 1 1 i r 1 1 1 1 1 i r 
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(b) V?/, = 0.10, Se = 0.05 foot, = 3.955 x 10 -5 square foot. 


Figure 12. - Effect of turbulence on mixing in first pure-standing 
tangential mode. Oscillation frequency, 1000 cycles per second; 
stream velocity, 500 f.eet per second; axial station, 1.25 feet; 
fraction .of cycle, 0.5. 
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Figure 13, - Mixing in first pure- standing tangential mode vith injection at points 
away from duct axis. 
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Combustor over- all equivalence ratio 
("based on unburned inlet air) 


Figure 14. - Typical combustor screech 
limits (unpublished data). 
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